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Lateral imaging of the superconducting vortex lattice using 
Doppler-modulated scanning tunneling microscopy 
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By spatially mapping the Doppler effect of an in-plane magnetic field on the quasiparticle tunneling spectrum, 
we have laterally imaged the vortex lattice in superconducting 2i^-NbSe2. Cryomagnetic scanning tunneling 
spectroscopy was performed at 300 mK on the a6-surface oriented parallel to the field H. Conductance images 
at zero bias show stripe patterns running along with the stripe separation varying as H~^-^. Regions of 
higher zero-bias conductance show lower gap-edge conductance, consistent with spectral redistribution by 
spatially-modulated superfluid momentum. Our results are interpreted in terms of the interaction between 
vortical and screening currents, and demonstrate a general method for probing subsurface vortices. 

PACS numbers: 74.55. +v, 74.25.Uv, 74.70.Ad 



In response to an applied magnetic field, type-II su- 
perconductors experience a diamagnetic current that cir- 
culates along the sample edge. Above the lower critical 
field, field can penetrate into the superconductor via a 
lattice of vortices, each consisting of a paramagnetic cur- 
rent loop enclosing a flux quantum.— The vortex lattice 
can be imaged using techniques sensitive to variations 
in the local magnetic field such as Bitter decoration^ 
and Lorentz microscopy,^ or with scanning tunneling mi- 
croscopy (STM) which probes the local quasiparticle den- 
sity of states (DOS). STM imaging of the vortices is pos- 
sible by virtue of bound states and suppressed supercon- 
ducting gap in the vortex corej^*^ Because of this reliance 
on vortex-core states, STM imaging has been largely lim- 
ited to the cross-sectional geometry, i.e. with the vor- 
tices piercing the sample surface. An earlier STM study 
versus field direction has shown the density of vortices 
to decrease as the field is tilted away from the surface 
normal^^ii For fields parallel to the surface, the vortex 
cores become buried in the bulk, making them difficult 
to probe directly. In this lateral field geometry, vortex 
lattices have been imaged by Lorentz microscopy, but 
only in highly 2D superconductors where pancake vor- 
tices decorate in-plane flux lines ^^i^ 

In this letter, we report on lateral imaging of the su- 
perconducting vortex lattice using cryomagnetic scan- 
ning tunneling spectroscopy. By mapping the zero-bias 
tunneling conductance over the a6-surface of supercon- 
ducting 2i^-NbSe2 in an in-plane magnetic field and at 
300 mK, we have observed distinct stripe patterns whose 
orientation and spacing versus the field can be directly 
attributed to the in-plane flux lattice. Our observa- 
tions are interpreted in terms of the interaction between 
the diamagnetic screening current and the paramagnetic 
vortical currents, which results in a spatially-modulated 
Doppler effect on the quasiparticle DOS spectrum. 

The STM used in our experiment was specially de- 
signed for the magnetic field to be applied parallel to 
the sample surface, as shown in Fig. [D^e). The STM 



is mounted inside a ^He cryostat which is inserted into a 
superconducting solenoid. The Pt-Ir tips used were field- 
emitted in situ to ensure stable vacuum tunneling, and 
RF-filters were used throughout the wiring to maximize 
the spectral resolution. The dl/dV conductance spec- 
tra were acquired by lock-in amplification with a 20 /iV 
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FIG. 1. Lateral imaging of the superconducting vortex lattice 
in 2i7-NbSe2 at various fields and 300 mK. Panels (a) to (c) 
show 380 X 380 nm^ spatial maps of the normalized zero-bias 
tunneling conductance Go for fields of 0.09 T, 0.25 T, and 
0.4 T, applied along the [100] direction (white arrow in (a)). 
Panel (d) shows the average Go along the direction perpen- 
dicular to the stripes in the 0.09 T data. Panel (e) shows a 
schematic of our experiment, with the STM measuring the 
a6-surface, across which flows a diamagnetic screening cur- 
rent (red line). In the heuristic model discussed in the text, 
paramagnetic currents (blue loops) circulating the subsurface 
vortices perturb the screening current, thus spatially modu- 
lating Go and producing the observed Go image contrast. 
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FIG. 2. Separation between the centers of the Go stripes, for 
fields between 0.09 T and 0.4 T, plotted on a log-log scale. 
The data (circles) are fitted to oc if-o-54±o.o9 ^^^^ ^-^^^^^ 
good agreement with the expected H~^'^ dependence of the 
vortex lattice parameter on field. 



excitation at 505 Hz, and the typical high-bias junction 
impedance was 10 MQ. 

Single crystals of 2i^-NbSe2 were grown by an iodine 
vapor transport technique. The crystals had critical 
temperatures of ~ 7.2 K, and upper critical fields of ~ 
5 T and 15 T respectively for field perpendicular and 
parallel to the a6-plane. The crystals measured were ^ 
5 X 5 X 0.5 mm^ in size, with the wide faces normal to 
the c-axis. The crystals were oriented by X-ray diffrac- 
tion, cleaved just before being loaded into the STM and 
cooled to 300 mK in zero field. STM topography revealed 
atomically smooth surfaces, with hexagonally arranged 
Se atoms modulated by triangular charge density waves. 
The magnetic field was applied along the [100] direction 
with ~ 2° precision. 

Figure [1] shows spatial maps of the zero-bias conduc- 
tance. Go- The Go data was normalized relative to the 
above-gap conductance at 4 mV. Panels (a) to (c) show 
plots of the data for 0.09 T, 0.25 T and 0.4 T respectively. 
Regions of high Go, hereafter referred to as stripes, can 
be seen in each plot, running parallel to the in-plane field 
and spaced at regular intervals. Crystals from different 
growth batches showed the same stripe patterns, thus at- 
testing to their general reproducibility. The average half- 
width of the stripes is 35±15 nm over the field range mea- 
sured. This is comparable to the expected size of a vortex 
core 2^ab ^ 20 nm, where is the zero-temperature su- 
perconducting coherence length in the a6-plane for 2H- 
NbSe2 The brightness of the stripes alternates between 
adjacent stripes, an effect more clearly seen in panel (d), 
which shows a profile plot of the average Go along the 
direction perpendicular to the stripes in the 0.09 T data. 

To interpret the stripe patterns as a manifestation of 
the subsurface flux vortices, we consider quantitatively 
how these patterns vary with the applied field H. Figure 
[2] shows a log-log plot of the separation between stripe 
centers over the field range measured in this experiment. 
The stripe separation decreases with increasing field and 
is fitted to oc i;f-0-54±o.09^ good agreement with the 




FIG. 3. (a) Series of normalized dl/dV spectra measured as 
the tip is scanned perpendicularly across a stripe at 0.1 T and 
300 mK. The midgap states rise over a ~ 35 nm wide region, 
thus accounting for the stripe patterns seen in the Go im- 
ages, (b) Spatially-averaged dl/dV spectra over the center of 
a stripe (red curve) and just between two stripes (blue curve). 
The zero-field spectrum (green curve) is included for compar- 
ison, to show the overall effect of the applied field. Over the 
stripe, the zero-bias dl/dV is higher while the coherence-peak 
dl/dV is lower, an effect clearly seen in the difference of the 
two curves, as shown in (c). 



expected H~^-^ dependence of the vortex lattice param- 
eter on field. It is important to note that the stripe sep- 
aration we observed is 0.86±0.10 times the Abrikosov 
lattice parameter. This observation can be qualitatively 
explained by identifying every other stripe with the row 
of fiux lines closest to the surface, and the in-between 
stripes as coming from the next closest row; this picture 
would be consistent with the alternating brightness be- 
tween adjacent stripes noted above. A more quantitative 
explanation of the stripe separation observed would need 
to consider vortex lattice distortions due to the super- 
conducting anisotropy of 2iJ-NbSe2 over the field range 
we measuredJ^ 

To visualize spectroscopically how the stripe patterns 
emerge, we analyze detailed variations in the dl/dV spec- 
trum across a typical Go image. Figure [3Ka) shows a se- 
ries of dl/dV spectra measured while the tip is scanned 
perpendicularly across a stripe region at 0.1 T. It is clear 
that the midgap states rise appreciably over a ~ 35 nm 
wide region, thus accounting for the stripe patterns seen 
in the Go image described above. Panel (b) shows the 
spatially- averaged spectra over the center of a stripe and 
just between two stripes. The zero-field spectrum is in- 
cluded for comparison, to show the overall effect of the 
applied field. Comparing the two curves at 0.1 T, the 
zero-bias dl/dV is higher while the dl/dV near the co- 
herence peak at ~ ±1.4 mV is lower, over a stripe than 
between stripes. This effect can be seen more clearly in 
the difference of the two spectra, as shown in panel (c). 
Here we emphasize the subtlety of this midgap-states ef- 
fect, i.e. the difference in Go is ~ 5 nS, requiring high 
measurement sensitivity to detect. 
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To understand the physical origin of the spectral evolu- 
tion versus field, we consider the recent theoretical model 
of Zhang et al. (Ref. [isl), which calculates the quasi- 
particle DOS spectrum for a current-carrying supercon- 
ductor. Essentially, the super current adds a Doppler 
term proportional to Vp ■ Qs to the quasiparticle en- 
ergy dispersion E^^ where Vp is the Fermi velocity and 
Qs is the superfluid momentum. This Doppler shift in 
/c-space causes a redistribution of the energy-gap size 
along the Fermi surface (FS), thus modifying the quasi- 
particle DOS spectrum. With increasing qr^, the co- 
herence peaks are suppressed in height and shifted to 
higher energies J^i^ When the Doppler term becomes 
sufficiently large relative to the gap amplitude, parts of 
the FS become gapless, causing the zero-energy DOS to 
increase with the field. All of these spectral behaviors, 
namely suppression of the coherence peak and enhance- 
ment of the midgap states, are observed in our tunneling 
measurements, thus directly implicating the Doppler ef- 
fect in the appearance of finite Gq with in-plane field. 

Finally, we consider the interaction between the vorti- 
cal and screening currents, in order to explain why the 
Doppler redistribution of the quasiparticle DOS spec- 
trum shows spatial modulation that is correlated with 
the vortex lattice. These two currents flow in oppo- 
site directions near the sample surface, as shown in 
Fig. me), and the counter-flowing current lines would 
be denser over certain regions in order to avoid crossing 
each other, thereby locally enhancing the supercurrent 
density and thus the superfluid momentum. A spatially- 
modulated superfluid momentum implies a spatially- 
modulated Doppler effect, which could then produce the 
stripe patterns seen in our in-field STM measurements. 
In this heuristic model, the stripes in the Go image corre- 
sponding to high-Doppler regions are over vortices, while 
the spaces between stripes corresponding to low-Doppler 
regions are between vortices. Although this model is 
physically plausible, it should be noted that there are 
alternative models indicating the opposite scenario, i.e. 
larger current density and thus stronger Doppler effect 
between vortices than over vortices ^^i^ Applicabilities 
of these models depend on quantitative details of the 
model assumptions, in particular the field strength rela- 
tive to the thermodynamic critical field and the coherence 
length relative to the London penetration depth. Further 
theoretical work is needed to more rigorously justify our 
model and quantitatively interpret our data. Here it is 
also worth noting that the subgap features at ^ ±0.7 mV 
seen in our spectra could be due to the multiband nature 
of 2i^-NbSe2, as manifestations of a smaller supercon- 
ducting energy gap which is also Doppler-redistributed 
by the in-plane field^i^— Elucidation of this scenario 
would require a multiband generalization of the Zhang 
et al. modelJ^, but would not affect the stripe patterns 



we observed or our interpretation of their physical origin. 

In summary, we have used Doppler-modulated scan- 
ning tunneling microscopy, performed at 300 mK, to lat- 
erally image the in-plane vortex lattice in superconduct- 
ing 2i^-NbSe2. This technique can in principle be applied 
on various other superconductors to probe the length 
scales and spatial symmetry of subsurface flux vortices. 
It may also potentially be used to study other fundamen- 
tal vortex phenomena, such as single-vortex entry and 
vortex-surface interactions. 

We acknowledge useful discussions with J. R. Clem, 
C.-R. Hu, V. G. Kogan and A. Paramekanti. This work 
was supported by NSERC, CFI-OIT and the Canadian 
Institute for Advanced Research. Part of this work was 
carried out at BNL, which is operated for the U.S. De- 
partment of Energy by Brookhaven Science Associates 
(Grant No. DE-Ac02-98CH10886). 

^A. A. Abrikosov, Sov. Phys. JETP 5, 1174 (1957). 

^U. Essmann and H. Truble, Phys. Lett. A 24, 526 (1967)] 

"^T. Matsuda, O. Kamimura, H. Kasai, K. Harada, 

T. Yoshida, T. Akashi, A. Tonomura, Y. Nakayama, 

J. Shimoyama, K. Kishio, T. Hanaguri, and K. Kitazawa, 

[Science 294, 2136 (2001), 

Caroli, P. G. D. Gennes, and J. Matricon, 

Phys. Lett. 9, 307 (1964)] 
^H. F. Hess, R. B. Robinson , R. C. Dynes , J. M. Valles, and J. V. 

Waszczak, Phys. Rev. Lett. 62, 214 (1989) 
^ H. F. He ss, G. A. Murray, and J. V. Waszczak, 
jPhys. Rev. Lett. 69, 2138 (1992)| 

'TL R Hess, C. 'K. Murray, and J. V. Waszczak, 

Phys. Rev. B 50, 16528 (1994)1, 
^V. K. Vlasko-Vlasov, A. Koshelev, U. Welp, G. W. Grabtree, 

and K. Kadowaki, Phys. Rev. B 66, 014523 (2002) 
^M. Yasugaki, K. Itaka, M. Tokunaga, N. Kameda, and 

T. Tamegai, Phys. Rev. B 65, 212502 (2002) 
^ "G. S . Oglesby, E. Bucher, G. Kloc, and H. Hohl, 

[Journa l of Grystal Growth 137, 289 (1994) 
^-' ^F. D. Gallaghan, M. Laulajainen, G. V. Kaiser, and J. E. Sonier, 
jPhys. Rev. Lett. 95, 197001 (2005) 

^ ^L. J. Gampbell,~ M. M. Doria, and V. G. Kogan, 

iPhys. Rev. B 38, 2439 (1988)[. 
^■ ^D. Zhang, G. S. Tin g, and G.-R. Hu, 

JPhys. Rev. B 70, 172508 (2004) | 

Fulde, m Tunneling Phenomena in Solids, edited by 

E. Burstein and S. Lundqvist (Plenum, New York, 1969) p. 427. 
i^A. Kohen, T. Proslier, T. Gren, Y. Noat, W. Sacks, H. Berger, 

and D. Roditchev, Phys. Rev. Lett. 97, 027001 (2006) 
i^E. H. Brandt, Journal of Low Temperature Physics 42, 557 (198T)] 
I'^A. D. Hernandez and A. Lopez, Phys. Rev. B 77, 144506 (2008) 
-•^^E. Boaknin, M. A. Tanatar, J. Paglione, D. Hawthorn, F. Ron- 

ning, R. W. Hill, M. Sutherland, L. Taillefer, J. Sonier, S. M. 

Hayden, and J. W. Brill, Phys. Rev. Lett. 90, 117003 (2003) 
i^G. L. Huang, J.-Y. Lin, Y. T. Ghang, G. P. Sun, H. Y. 

Shen, G. G. Ghou, H. Berger, T. K. Lee, and H. D. Yang, 

Phys. Rev. B 76, 212504 (2007) 
^' ^J. G. Rodrigo and ^S. Vieira, 

[Physica G: Supercond uctivity 404, 306 (2004 )^ 
^ ^L Guillamon, H. Suderow, F. Guinea, and S. Vieira, 

|Phys. Rev. B 77, 134505 (20087) 



